Breakdown of Hydrodynamic Transport Theory in the Ordered Phase of Helimagnets 
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It is shown that strong fluctuations preclude a hydrodynamic description of transport phenomena 
in helimagnets, such as MnSi, at T > 0. This breakdown of hydrodynamics is analogous to the 
one in chiral liquid crystals. Mode-mode coupling effects lead to infinite renormalizations of various 
transport coefficients, and the actual macroscopic description is nonlocal. At T — these effects 
are weakened due to the fluctuation-dissipation theorem, and the renormalizations remain finite. 
Observable consequences of these results, as manifested in the neutron scattering cross-section, are 
discussed. 

PACS numbers: 



Soft modes, or excitations at low frequencies and small 
wave numbers, in many-body systems are crucial for the 
behavior at long distances and times. They also lead to 
characteristic nonanalyticities in the temperature (T) de- 
pendence of observables that are given by integrals over 
the excitation spectrum, such as the specific heat. Exam- 
ples include shear diffusion in classical fluids Q], the 'dif- 
fuson' in quenched disordered electron systems at T = 
, magnons in magnetically ordered solids Q , the direc- 
tor fluctuations in liquid crystals 0], and many others. 
Interactions between such soft modes, often referred to 
as mode-mode coupling effects, lead to even more inter- 
esting phenomena. In classical fluids, they produce the 
so-called long-time tails, or non-exponential decay of time 
correlation functions, e.g., the transverse velocity auto- 
correlation function [5j . In smectic and cholesteric liquid 
crystals they are even stronger, and lead to a breakdown 
of hydrodynamics, with certain viscosities diverging for 
vanishing frequency w as 1/w 6]. The weak-localization 
anomalies in disordered electron systems Q belong to the 
same category of effects, although they are usually not 
cast in this language. Technically, the mode-mode cou- 
pling effects take the form of nonlinearities in dynamical 
equations, or loops in a field-theoretic description. 

The above examples make it clear that mode-mode 
coupling effects, and the singularities induced by them, 
are relevant for both classical and quantum many-body 
systems. An obvious question is therefore what happens 
to the classical singularities as T — > 0. From the structure 
of the classical singularities it is likely that they become 
weaker with decreasing temperature, while quantum sin- 
gularities may develop that become stronger, but this has 
never been investigated. 

One purpose of the present Letter is to study an ex- 
plicit example of what happens to classical mode-mode 
coupling singularities as T — > 0. Our example are 
helimagnets, such as MnSi, which have recently been 
shown to exhibit interesting analogies with smectic and 
cholesteric liquid crystals. These analogies exist because 



the order parameters in the two systems are essentially 
identical and the symmetries of the respective coarse- 
grained free-energy functionals are closely related. Con- 
sequently, a tree-level analysis of both a phenomenolog- 
ical and a microscopic theory of helimagnets Q reveals 
a Goldstone mode, and its manifestation in observables, 
that are very similar to the corresponding properties of 
chiral liquid crystals 0, @- These properties are very 
unusual because of the anisotropic nature of the helical 
Goldstone mode, which is softer in the direction trans- 
verse to the pitch vector of the helix than in the longitu- 
dinal direction. 

A second purpose of this Letter is to investigate 
whether these analogies extend beyond the tree level, i.e., 
whether the mode-mode coupling effects in helimagnets 
are as strong as in chiral liquid crystals. We find that 
this is indeed the case. Some of the most interesting 
and dramatic effects are the fluctuation renormalizations 
of the transport coefficients that determine, for exam- 
ple, the line width of the neutron scattering signal due 
to the helical Goldstone mode. At T > we find that 
fluctuation effects are so strong that they lead to infinite 
renormalizations of the transport coefficients and elastic 
constants. Physically, this implies that at long length 
scales, and for low frequencies, the actual macroscopic 
description is nonlocal. In the quantum or zero tempera- 
ture limit we find that the divergencies are removed due 
to the structure of the fluctuation-dissipation theorem, 
and the transport coefficients and elastic constants sat- 
urate, albeit at a very low temperature. As we discuss 
below, our predictions will be observable in precision neu- 
tron scattering experiments. 

Our starting point is a phenomenological hydrody- 
namic theory [fj, for fluctuations in helimagnets. 
We will take into account the most important nonlinear 
terms, and will show that the nonlinearities qualitatively 
affect the physics. Let the mean magnetization (m) lie 
in the x-y plane, with the axis of the helix pointing in 
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the z-direction: 



(m)(x) = mo (cos qz, singz, 0), 



(1) 



with q the pitch wave number of the helix. Typically, 
q <C kp, with k-p the Fermi wave number. The fluctua- 
tions about this average state are described in terms of 
a dimensionless generalized phase variable u(x), which 
describes the fluctuations of the x and y-components of 
m, and the z-component of the magnetization, m z (x) = 
m(x) 0- The course-grained Hamiltonian for these vari- 
ables takes the form 



H[u,m] = — J dx [aq 2 m 2 — 2hm\ 



dx 
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Here a is the coefficient of the square-gradient term in 
a purely magnetic Hamiltonian, h is the magnetic field 
conjugate to m, c z and c± are elastic constants, and 
= dl + d 2 . The first line in Eq. J2J is just the Gaus- 
sian part of a Landau Hamiltonian for m; higher powers 
of m are not important for our purposes. The second line 
is the effective Hamiltonian for u, which is constrained by 
symmetry considerations 00. Translational invariance 
implies that only gradients of u can enter, while rota- 
tional invariance requires that pure rotations of the helix 
do no cost any energy, which restricts the combinations 
of gradients that can enter the Hamiltonian. Up to quar- 
tic terms in u, Eq. J2J is the most general Hamiltonian 
consistent with these constraints. 

Equation J2J completely describes the thermodynam- 
ics of the system via the partition function Z = 
Tr exp(— H/T) (we use units such that /cb = !)• The 
dynamics are described by coupled dynamical equations 
that follow 7] from the usual time-dependent Ginzburg- 
Landau theory for the magnetization 9] and read 
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Here 7 is a constant, and so are the transport coefficients 
T' u z x 11] . h is the (magnetic) field conjugate to to, and 

and C are Langevin forces. In Eq. i|3"aj) the dissipative 
term oc T^, as well as the Langevin force turn out 
to be unimportant for the long-wavelength effects con- 
sidered here, and we will neglect them in what follows. 

Consider the magnetic susceptibility or linear response 
function \i defined by the response of the average mag- 
netization to an applied magnetic field, 



(m(k,uj)) = x{k,Lo) h(k,w) 



(4) 



to the imaginary part or spectrum of x, x"(k, us), by the 
fluctuation-dissipation theorem [l^ | 



F mm {k,u) = X "{k,uj) coth(hw/2T). 



(5) 



The renormalizations of the elastic coefficients and trans- 
port coefficients are determined by examining how fluc- 
tuations modify the linear response function. To zeroth 
order in the non-Gaussian terms in Eq. J2J the latter is 
easily found, from Eqs. to be 



fll(k) - iLuT (k) 



aq 2 [-uj 2 + nl(k) - iLuT {k)}' 



Here 



fio(fc) =7a 1/2 9 [czk 2 z + cjfci] 



1/2 



is the hclimagnon resonance frequency, and 

r (fc) = iy^ + r x fci, 



(6a) 



(6b) 



(6c) 



with Y z .j_ = aq 2 T' z j_, is a damping coefficient. The spec- 
trum of x is proportional to the neutron scattering cross- 
section, and the renormalized counterpart of Tq deter- 
mines the line width of the helimagnon resonance |l3l |. 
The (tin) and (mu) correlation functions can be related 



(7) 



to F mm by Fourier transforming Eq. (|3a|) to obtain 
u(k,uj) — 7 [aq 2 m(k,uj) — h] /(—ioj). 



Loop corrections to xo due to the non-Gaussian terms 
in Eq. @, which lead to nonlinearities in the dynamical 
equations, can be calculated by means of standard per- 
turbation theory for x- The two one-loop diagrams are 
shown schematically in Fig.^ It is convenient to consider 





FIG. 1: One-loop contributions to the response function \- 
The two diagrams result from the cubic and quartic in u 
terms, respectively, in Eq. ©• 

the renormalization of Xq 1 , whose one-loop correction we 
denote by Sx^ 1 ■ We obtain 



where 



to 



[n§(*0- iwr (fc)]' 



[7 3 (fe,cj)+J 4 (fe,a;) 
-C(fe)] , (8) 



C(k) = TK 2 k\ 



The Langevin force Q is fixed by requiring that the cor- 
relation function F mm {k,bj) — (\m(k,uj)\) [12( is related and 
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Here 



transport coefficients for small T± are 



W 3 (fc,p) 
Wk(fc,p) 



[^3(*!,p) + ^3(fc,fc-p)] 

x[V 3 (p,k) + V 3 (p,p-k)}, (11a) 
V4(fc,fc,p) + 14(fe,p,fe) 
+V i (k,-p,p), (lib) 



with 



V 3 (fc,p) 
V 4 (fc,pi,p 2 ) 



(k z p± + 2p z k ± ) ■ (fei 
[{k-pi±)p2± + (pi± 
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The first and second terms on the right-hand side of Eq. 
iJHJ correspond to the two diagrams shown in Fig.^ The 
third term arises from a counterterm that needs to be 
added to the Hamiltonian in order to ensure that one 
expands about a proper saddle point, rather than its 
mean- field approximation |l4j |. A is an ultraviolet cut- 
off common to all momentum integrals. The technical 
effect of the counterterm is to subtract the k\ contri- 
butions that arise in perturbation theory and to ensure 
that the resonance frequency retains the structure given 
in Eq. (|6bl) . which is dictated by rotational symmetry. 

Using Eqs. iloKill in Eqs. ((H) - fTT|) leads to the leading 
order fluctuation renormalization of the response func- 
tion. Note the difference between the classical and quan- 
tum cases: In the classical limit, the coth in Eq. (JSJ 
contributes a factor T/{u> — lo\) to the integrand in Eq. 
IjlOajl . while in the quantum limit it contributes a step 
function. The important implication is that in the quan- 
tum limit, the analog of classical mode coupling effects 
are relatively less singular at low frequencies. 

Equation (JHJ) effectively defines a renormalization of 
the elastic coefficients and the transport coefficients [lij . 
Denoting the corrections to fl^ and To by 5fl 2 and 5T, 
respectively, we find 



sn 2 (k) 

5T(k,uj) 



1 
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aq z 



Re[/ 3 (fe,0)+7 4 (fe,0)-C , (fe)], (12a) 
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(12b) 



Equations l|12|l are the central results of this paper. In 
the classical limit, ft = 0, at zero external wave number 
and finite frequency the perturbative corrections to the 



5T Z = 



ST i 



128aq 4 c/ 2 M 
{^aq 2 ) 2 cj 2 
32iraq 4 c 1 / 2 T± 



T ln(A 2 r x /M) 



(13a) 
(13b) 



For the one-loop corrections to the elastic coefficients at 
zero frequency we find 



5c z 



5c± = 



3/2^ 
Cz T 

32nq 2 c 3 / 2 
cl /2 T 



ln(A 2 /fci) 
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(14a) 
(14b) 



Equations l|14|) can also be obtained from equilibrium sta- 
tistical mechanics considerations, by renormalizing the 
Hamiltonian, Eq. Note that T z , T±, and c± are 

renormalized to infinity, while c z is renormalized to zero 
[lif . Also note that the fluctuation corrections are all lin- 
early dependent on the temperature. This is due to the 
fact that all of the excitations are bosonic in character. 

Physically, the Eqs. (|14fl imply that the helimagnon 
state is, strictly speaking, not stable at finite tempera- 
tures. The same is true for certain liquid crystals that 
show one-dimensional order, and it shows that the analo- 
gies between helimagnets and liquid crystals pointed out 
m Ref. at a Gaussian level do indeed carry over to the 
respective nonlinear theories. This is a nontrivial obser- 
vation, since the dynamical equations are different in the 
two cases. However, practically speaking one must go to 
enormous length or time scales for these divergencies to 
dominate the zero-loop effects. In liquid crystals the sin- 
gularity is typically cut off by finite system size effects; 
in a helimagnet, it will be cut off by crystal-field effects 
that break the rotational symmetry. The strongest ef- 
fect is the renormalization of T z , which shows a linear 
divergency rather than a logarithmic one. 

We now turn to the quantum or zero-temperature 
limit, where the coth in Eq. JSJ turns in to a step func- 
tion. ST±, 5c Zl and 5c± are then finite by power counting. 
5T Z one might expect to still be logarithmically divergent, 
but a detailed inspection of the relevant integral shows 



4 



that it, too, is finite. We find 



sr z 



(jaq 2 ) 2 c 3 J 2 



(T = 0) 



(15) 



512 aq A cJ 

Comparing with Eq. (|14af> we see that 8T Z saturates at 
a value corresponding to T = |w|/4. All of the classical 
divergencies thus vanish at T = 0. 

Let us finally provide order-of-magnitude estimates for 
the magnitude of these fluctuation effects. The coeffi- 
cients a (which dimensionally is a length squared) and j 2 
(which dimensionally is an energy times a length cubed) 
have nothing fundamental to do with magnetism and are 
of order a ~ l/fc F and j 2 ~ £fAf, respectively. We de- 
note the Fermi energy, wave number, and velocity by cf, 
fcp, and vf, respectively, and the Stoner gap or exchange 
splitting by A. The tilde indicates that we are leaving 
out a numerical prefactor that one expects to be close 
to unity. For the elastic coefficients one expects c z ~ 
c±q 2 oc A 2 /wp. However, the model calculation of Ref. 7 
yielded a small prefactor on the order of 10 -3 for these 
coefficients, and we will use c z = 2c±q 2 ~ 0.001 X 2 /vf 
for our estimates. The transport coefficients V ' z and 
have again nothing to do with magnetism, and we expect 
~ T 1 ^ ~ VpT, with t the electronic elastic mean- free 
time. This leads to T z ~ Tj_ ~ (q/k F ) 2 VpT. 

Now consider the singularities at T > 0. In a solid- 
state system they are cut off at a length scale L that is 
inversely proportional to the spin-orbit coupling squared, 
or, since q is proportional to the spin-orbit coupling, 
L ~ 2-Kk^jq 2 [tJ. The cutoff frequency for the dynam- 
ical singularities is given by Qo(k 
estimates we obtain 



1/L). With these 



r z 



T e F L 
T 



0.003 • 



(16a) 



ln(fcpL), (16b) 



with i = v-pT the elastic mean-free path. For the elastic 
constants we find 

^~^i~-100^(f) 2 f Mk F L), (17) 

Parameter values for MnSi are (Q and references therein) 
e F « 23, 000 K, q/k F w 0.024,' and L « 7, 500 A. The 
value of cf/A is less certain, we will assume £f/A = 4, 
which leads to 6c z /c z « 0.1 at T — 10K. The correc- 
tions to the damping coefficients are much smaller. With 
a mean-free path I w 5, 000 A [Tjj one finds 8T Z /T Z w 
10 5 (5r ± /r ± w 10" 3 at T = 10 K. In less clean samples 
the fluctuations of T z and T± will obviously be larger. 

The temperature scale where the divergent classical 
results cross over to the finite expressions in the quantum 
limit follows from Eqs. I|13all5l) . and is corroborated by 
an inspection of the relevant integrals. It is given by 



With the same parameter values as above we find T* f» 
0.3 mK. In interpreting all of these quantitative results 
one needs to keep in mind that these are very rough esti- 
mates that depend strongly on some of the parameters. 

In conclusion, we find anomalously large fluctuation 
effects in the ordered phase of helimagnets which lead to 
a breakdown of hydrodynamics. In MnSi, the largest ob- 
servable signature is predicted to be a linear-in-T shift in 
the resonance frequency squared, with an absolute value 
of the shift of about 10% at T = 10 K. This should be 
observable by neutron scattering. 
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